N79-20036 


LOW SPEED AIRFOIL DESIGN AND ANALYSIS 6 

Richard Eppler 

Universitat Stuttgart, West Germany 
Dan M. Somers 

NASA Langley Research Center 


SUMMARY 

A low speed airfoil design and analysis program has been developed which 
contains several unique features. In the design mode, the velocity distribu- 
tion is not specified for one but many different angles of attack. Several 
iteration options are included which allow the trailing edge angle to be 
specified while other parameters are iterated. For airfoil analysis, a panel 
method is available which uses third-order panels having parabolic vorticity 
distributions. The flow condition is satisfied at the end points of the 
panels. Both sharp and blunt trailing edges can be analyzed. The integral 
boundary layer method with its laminar separation bubble analog, empirical 
transition criterion, and precise turbulent boundary layer equations compares 
very favorably with other methods, both integral and finite-difference. 
Comparisons with experiment for several airfoils over a very wide Reynolds 
number range were very favorable. Applications to high lift airfoil design 
were also demonstrated. 


INTRODUCTION 

The application of potential flow theory together with boundary layer 
theory to airfoil design and analysis was accomplished many years ago. Since 
that time, high speed computers have allowed results to be obtained more 
cheaply and quickly than through the use of wind tunnels. Accordingly, the 
tendency today is toward more and more conwonly applicable computer programs. 
The programs reduce the amount of required wind tunnel testing to that of 
fundamental phenomena and allow airfoils to be tailored to each specific 
application. 

The program described in this paper has been developed over the past 20 
years. We hope to demonstrate that it has reached a stage where some progress 
can be made in low speed airfoii design and analysis. 

This paper does not repeat all of the details included in other papers. 
Special emphasis is given, however, to those features which are new or 
different from those in other matnematical models. 
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skin-friction coefficient 
airfoil chord 

section profile-drag coefficient 
section lift coefficient 

section pitching-moment coefficient about quarter-chord point 
shape factor, 

^2 

shape factor, -- 

62 


length 
poi nt 

Reynolds number based on free-stream conditions and airfoil chord 

number based on local conditions and boundary layer 
momentum thickness ^ 

thickness of airfoil 

free-stream velocity 

local velocity on airfoil 

velocity; main pressure recovery design variable 

airfoil abscissa; length 


angle of attack relative to zero-lift line, deg 

angle of attack relative to zero-lift line for veloci 
in design method, deg 

boundary layer displacement thickness 

boundary layer momentum thickness 


ty specification 


boundary layer energy thickness 
kinematic viscosity 


total amount of pressure recovery 
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u 


\ 




initial slope of pressure recovery 
Subscripts: 

^ lower surface; local point on airfoil 

im lower surface main pressure recovery 

N last point on airfoil 

n variable number 

u upper surface 

urn upper surface main pressure recovery 

« free-stream conditions 


DISCUSSION 
Inviscid Method 


The potential flow part of the mathematical model is incompressible at 
this time. Two different modes of operation are available. 

Design mode .- The first mode is the inverse or design method described 
'•n references 1 and 2. This method differs from other inverse methods in that 
the velocity distribution is not specified for only one angle of attack. 
Instead, angles of attack which will result in constant velocity over speci- 
fied segments of the airfoil are input. In other words, pairs of parameters 
are specified: the first being the segment of the airfoil; the second, the 
angle of attack relative to the zero-lift line, a*, which will result in 
constant velocity over that segment. (See fig. 1.) Of course, some matching 
conditions must be met to guarantee a smooth velocity distribution for all 
angles of attack. Toward the trailing edge, on both surfaces, a main pressure 
recovery can be specified. Finally, a short closure contribution must be 
introduced to insure that the trailing edge will be closed. The example air- 
foil shown in figure 1 is specified by the following: 

(1) For the uoper surface segment from the trailing edge forward to about 
x/c = 0.15, a* = 8°. Within this segment, the main pressure recovery is 
specified starting at x/c = 0.50. 

(2) For the upper surface segment from x/c = 0.15 to 0.05, a* = 10°. 

(3) For the upper surface segment from x/c = 0.05 to the leading edge, 

a* = 12°. 

(4) For the entire lower surface, a* = 2°. The main pressure recovery 

on the lower surface is specified starting at about x/c = 0.45. 
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In reality, the segnsnts corresponding to the various o* values are 
not specified In x/c but rather In the conformal mapping plane In which the 
airfoil Is j*epresented by a circle. So far, no difficulties have arisen In 
correlating the arcs of the circle with the segments of the airfoil. 

It should be remembered that for any given velocity distribution there 
does not necessarily exist a "normar airfoil. For example, the closure con- 
tributions could be quite large which would result In a very large trailing 
edge angle. The closure contributions could also give rise to a region of 
negative thickness near the trailing edge. Accordingly, several Iteration 
options have been Included which allow the trailing edge angle to be specified 
while certain o* valiies or the total amount of pressure recovery 1« Iterated. 
The choice of Iteration option allows questions such as the following to be 
answered: What laminar bucket width Is possible given a certain amount of 
pressure recovery? What amount of pressure recovery Is required to produce 
the desired laminar bucket width? The Iteration option selected for the 
example airfoil Iterates the amount of upper and lower surface pressure recov- 
ery while holding the a* values fixed. 

Analysis mode. - The second mode of operation Is an airfoil analysis 
Tiethod. The method employs panels with distributed surface singularities 
(fig. 2). The panels are defined by a third-order spline fit of the airfoil 
coordinates with the end points of the panels being the Input airfoil coordi- 
nates themselves. The surface singularities consist of a parabolic vorticit.y 
distribution. The flow condition, which requires the Inner tangential 
velocity to be zero. Is satisfied at each airfoil coordinate (I.e., at the 
end points of the panels, not the mid-points). Thus, no restrictions are 
placed on the point distribution, no smoothing or rearranging of the coordi- 
nates Is performed, only the original airfoil coordinates are used. An 
option Is Included, however, by which additional points can be splined In 
between the original coordinates. This option Is helpful If a portion of the 
airfoil has a sparse number of points or If part of the airfoil Is to be 
geometrically rotated about a flap hinge point. In the latter cast, the 
connection between the forward portion of the airfoil and the flap Is defined 
by an arc consisting of additional points which are generated automatically 
according to an Input length. 

As In other panel methods, a singularity arises from the circulation 
around the airfoil which Is unconstrained unless required to satisfy a Kutta 
condition. Two different cases are Involved. 

The first case Is a sharp trailing edge having either a zero or non-zero 
angle. The Inner tangential flow condition falls In this case at the trail- 
ing edge. Therefore, It Is replaced by the condition that the normal 
velocity, relative to the bisector of the trailing edge angle, be zero. This 
condition can only be satisfied If the vorticitles on the upper and lower 
surfaces approach the same value but with opposite rigns as the trailing edge 
Is approached. This means that the velocities have to be the same on both 
sides of the trailing edge. Thus the normal flow condition along with equal 
velocities on both sides of the trailing edge can be considered as a 
Kutta condition. 
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Unfortunately, a second singularity is caused by a finite change in the 
airfoil which results in an infinite change in the 
velocity at the sharp trailing edge. No wonder the equation system is sinou- 
nnl 1°*^ 6y®'^y ^®tta condition. Thus, an additional equation is required while 

eouation^col!«;iitrif equation system can be omitted. This additional 

extrapolation of the vorticity to the trailing edge 
in'^tho at the trailing edge. The omission of one of the equations 

gSve?ned condition is 

Lactfv^iiaJlar hJp equation system is not 

s^'?9Ular, due to small numerical approximations. So, none of the 

quations IS omitted and the entire equation system is multiplied by the 

ly'as'pofsiMel^a are siJJed'L'ac^urate- 

very Precise! least-squares sense. The results of this procedure are 

The second case is a blunt trailing edge (fig. 3) - one in which the 

twS^HifJp surface trailing edge points are not the same. For this case 

two different procedures have been examined. ’ 

The first procedure extends the airfoil along straight “wake limits” 

opposite vorticity (fig. 3, top). At both the upper and lower 
points, a normal velocity condition must be satisfied 
No flow condition is satisfied in the wake. 

WOV.+/-! second procedure introduces at the base of the airfoil not only a 
vorticity but also a source distribution (fig. 3, bottom). Both are linear 

ititP^^H^rnr^ over the base length and determined such that no flow singular- 
itites occur at either of the two trailing edge points. 

Both blunt trailing edge models have one more flow condition than 
unknown vorticities. The n + 1 equations for n unknowns are treated in 

rHi-ira? trailing edge case even though they are much less 

critical with respect to circulation changes. 

For all cases the computing times are moderate; the results, very precise 
shlZTnfTgull I example ai-rfoil fr 


Viscous Method 

c-imnil^® boundary layer development is computed by a 

pnSfnnc^^^Tt “Sing, like many others, integral momentum and energy 

equations. It has been shown that laminar boundary layer development is pre- 
dicted quite well by this method. The turbulent boundary layer ?o^?ines Ire 

factor'll awe en^PT^ical skin friction, dissipation, and shape 

further errors are introduced by mathematical simplifications 
.Ike integrating the ordinary differential equations from the momentum and 
energy laws by averaging the right sides of the equations. 
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Of special interest are the predictions of transition from laminar to 
turbulent boundary layer and the separation of the turbulent boundary layer. 
The tendency toward separation is determined solely by the shape factor 



where 6^ is energy thickness and eg is momentum thickness. (Note that ^32 

has the opposite tendency from H^g which contains the displacement thick- - 

ness 6.J instead of the energy thickness.) For laminar boundary layers 

there exists a constant and reliable lower limit of H^g, which equals 1.515 

and corresponds to laminar separation. For turbulent boundary layers no such 
unique and reliable limit has been determined. It can be stated, however, 
that the turbulent boundary layer will separate if H^g goes below 1.46 and 

will not separate if H^g remains above 1.58. It has been noticed that 

thicker boundary layers tend to separate at lower H^g values. In the 

present method, turbulent separation is predicted if H^g drops to 1.46. 

This is a fairly good assumption because the method usually predicts relative- 
ly low values of H^g. The uncertainty is not as bad as it first appears in 

that Hjg changes rapidly near separation. Nevertheless, results must be 

checked carefully with respect to turbulent separation. 

The second feature of special interest is the prediction of transition. 

Two different procedures are in vogue today - the amplification method and 
the shape-factor —• Reynolds-number method. The first procedure require? much 
more computing time because many frequencies must be traced to find the one 
wave which is suddenly amplified to the ratio set as the transition limit. 

The better procedure cannot be selected until we know more about transition. 

The differences between the two procedures are small for normal airfoil appli- 
cations as the local Reynolds number changes quickly near transition. So, we 
still use the simple criterion shown in figure 5. Thus, the transition 
Reynolds number depends only on the shape factor H 22 * Adverse pressure 

gradients and, hence, low values of H^g result in lower transition Reynolds 

numbers and vice versa. As will be seen later, it is very informative to plot 
the boundary layer development in this form. 


It must be pointed out that the boundary layer development inmediately 
after transition has a significant influence on the entire flow. In our 
program the prediction of transi :ion results in a switch from the laminar 
skin friction, dissipation, and shape factor laws to the turbulent ones, 
without changing and 6-. This is also done if laminar separation is 


predicted before the transition criterion is reached. The H-g development 
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for the linearly decreasing velocity distribution defined by 
V = VJl . x/i) 

is shown in figure 6 for Reynolds numbers ranging from 0.125 x 10® to 32 x 10® 

^ illustrates the well-known fact 

that the laminar boundary layer shape factor H32 and laminar separation are 

transitioS^JLJc^h^J^ also shows that for high Reynolds numbers, 

transUion occurs befP.e laminar separation. £or turbulervM>ound»ry layers, 

H32 and separatiDff do depend on Reynolds number. The most important infor- 
mation to be gained from figure 6 is the behavior of H32 at the beginning 
of the turbulent boundary layer. For higher Reynolds numbers, H32 increases 

1.7. For lower Reynolds numbers this 
increase is less rapid and the max-"num values of H32 are lower. For a 

Reynolds number of 0.126 X 10«, remains below 1.58 which means that the 

rethod cannot determine whether or not an attached turbulent boundary layer 
exists. Such results must be studied in more detail. In figure 7 the laminar 
and turbulent skin-friction laws, C^(R Hj^) are presented™ The laS? 

?" "e™®"* The turbiTent law Is a slightly modified Ludwleq- 

tellSS r '®“ experimentally derived and 
ested for R^ between 10 and 105 as shown by the phantom lines In figure 7. 

Below R^^ equal to 10 , these lines continue in some manner. The flat plate 

fi!5 investigated in more detail and the results Indicate that the 

flat plate line continues more or less steadily and finally bends down to the 
_aminar law line rather steeply, depending on roughness and free-stream 
turDuience. 

In our method, the Ludwieg-Ti liman law is extrapolated along straight 
lines. This probably represents an upper limit for C. for R, less than 10^ 

But it is obvious that for R^^ equal to 10^ tht laminar and turbulent laws 

differ little and for lower R the turbulent C. values are below the 
laminar ones. 2 ' 


results such as vhose shown in figure 6, it is interesting to look at 
tne values computed by our method. For every point x/i. and H 

32 


62 

with R, 


is shown 


are known and, hence, is known. The variation of C. 

I 

6^5® result is remarkable. The curves for R = 0.125 x 10 ® and 
m u ? "2^ even come close to an area where one can confidently 
speak of a turbulent skin-friction law. Such "underdeveloped" turbulent 
boundary layers exist -.n nature only in the form of laminar separation bubbles 
After examining many different cases, it was determ® ed that there is a certain 
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analogy between the predicted ''underdeveloped'' turbulent boundary layers with 
low H 32 values and laminar separation bubbles. The boundary layer results 

computed by our method show a stronger bubble analogy as Reynolds number 
decreases and as the adverse pressure gradient after transition becomes 
steeper. If the analogy occurs in the results, the only way to alleviate it 
is to reduce the adverse pressure gradient after transition. The experimental 
results for laminar separation bubbles show the same tendencies. Accordingly, 
it 's very helpful to have this bubble analog in the computed results. 

In summary, th' 6 * b 0 Xindary layer method has generated good" results formany, 
very different cases. It should be noted, however, that no fundamental 
problem exists in replacing the boundary layer subroutines in our program with 
other subroutines. Some applications will be discussed in the next section. 

But before that a few comparisons with other boundary layer mechods will be 
made. 


Two shape factor developments computed by a program written by Konhauser, 
which uses the Cebeci-Smith method (ref. 4), are shown by dashed curves in 
figure 6 . The two curves, which are for the Reynolds numbers of 8 x 10® and 
16 X 10 ®, agree quite well with the present method for x/i up to about 0 . 2 . 

Then as x/i increases, the results computed by the Cabeci -Smith method show 
considerably less tendency toward separation. This demonstrates that the 
separation limit of equal to 1.46 that we use is conservative with 

respect to turbulent separation. Comparisons for lower Reynolds numbers are 
not possible at this time because Mr. Konhauser has been unable to obtain 
results from the Cebeci-Smith method at lower Reynolds numbers. 

The development of displacement thickness along the upper surface of an 
RAE 101 airfoil for an angle of attack of 8 .^® and a Reynolds number of 
1.6 X 10® is shown in figure 8 . This case was computed by J. L. Hess using a 
Cebeci-Smith program with different numbers of elements and different smoothing 
procedures (ref. 5). The results from the program as indicated by the symbols 
agree very well with those computed by Hess with the greatest number of 
elements. This agreement is remarkable indeed knowing that this is the first 
comparison of this type which we have made and Is not the result of a careful 
choice of data. 

The curves in figure 9 demonstrate that, in using integral momentum and 
energy laws, the introduction of mathematical simplifications can cause much 
larger errors than those which result from the use of a one parameter method. 

In figure 9, our method and a method developed by L. Truckenbrodt (ref. 6 ) are 
compared with experimental results obtained by Wortmann (ref. 7). Truckenbrodt 's 
method is based upon the same skin-friction and dissipation laws as the present 
method, but includes further mathematical simplifications which produce an 
error of about 50 percent for the adverse pressure gradient shown. 

in the present method, the momentum thickness at the trailing edge is 
used for the calculation of the drag by a Squire-Young type formula. We have 
found that our method predicts slightly higher drag values than those measured 
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experimentally. We hesitate to change the skin-friction laws or other parts of 
the program, however, as the differences between our predictions and experi- 
mental measurements depend upon the wind tunnel in which the experiments were 
performed. (See figs. 10 and 11.) 

It should be mentioned that the program includes a lift coefficient 
correction due to boundary layer separation but, as yet, does not include one 
due to boundary layer displacement thickness. 

As a final remark, the development of the laminar boundary layer should 
be discussed. It is very informative to plot this development as shown in 
figure 5. This plot reveals several important points. For a constant velocity 
segment (a = a*), the boundary layer approaches the Blasius solution having a 
shape factor of 1.573 and increasing momentum thickness 62. This 

corresponds to the vertical lines in figure 5. As the angle of attack is 
increased, the velocity distributions become concave over the fo»ward portion 
of the airfoil or, in other words, the airfoil "pulls a peak" at the leading 
edge. These concave distributions are similar to those which produce Hartree 
boundary layers. But, whereas the Haitree boundary layers result in lower but 
still constant shape factors, the curves in figure 5 show increasing H-,.. with 

increasing R, for a greater than a . This means that these velocity 
“2 

distributions are more concave than the Hartree (power law) distributions. 

Thus, as the angle of attack is increaseu even more, these distributions will 
result in laminar separation at the leading edge. This problem is eliminated 
by the introductioi; of segments having higher a* values near the leading 
edge. Obviously, it is much easier to control the level opment of the shape 
factor by manipulating o* values than by changing a given velocity distribu- 
tion at only one angle of attack. 

Appli rations 

In this section, we shall apply the mathematical model to a variety of 
airfoil problems. The final result is always a plot which includes c^ 

versus c^, c^ versus o, Cj^ versus o, and transition and separation 

versus c^ as is normally plotted for wind tunnel results. It is, of course, 

very easy to obtain more details such as pressure distributions, boundary 
layer development, and laminar separatio*’ bubble analogs. 

The first application is a sailplane airfoil designed for low drag at a 
Reynolds number of about 3 x 10° and a soft stall at a Reynolds number of about 
1 X 10°. The soft stall can easily be achieved by introducing a moderate 
concave pressure recovery on the upper surface and by preventing laminar 
separation and the rapid forward movement of transition with increasing angle 
of attack. The latter feature requires only increasing c* values toward 
the leading edge. 
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The theoretical results agree well with the experimental measurements 
obtained by D. Althaus as shown in figure 10. In the wind tunnel experiment? 
transition location was determined by the stethoscope method which seems to 
detect only fully developed turbulence and, thus, the experimental transition 
locations lie somewhat downstream of the theoretical ones although the trends 
with angle of attack agree well. The stall observed in free flight was very 
soft. 

The next application shows that the program produces reasonable results 
for higher Reynolds numbers as well. The coordinates of an NACA 642-618 

airfoil were input and the theoretical results are compared with the experi- 
mental measurements (ref. 8) in figure 11. 

The program can also be applied at very low Reynolds numbers. Airfoil 387 
was designed for model airplanes. At these low Reynolds numbers, the bubble 
analog indicated that only very slight adverse pressure gradients were poss.ble 
and, accordingly, a relatively thin airfoil (t/c = 0.09) resulted. This air- 
foil was recently tested by Volkers (ref. 9). The theoretical results compare 
favorably with experiment for a Reynolds number of 2 x 10® (fig. 12), even 
though the measurements do show the typical effect of laminar separation 
bubbles. For a Reynolds number of 1 x 105, the experiment shows even more the 
effect of laminar separation bubbles, but still with attached turbulent flow 
at the trailing edge. For a Reynolds number of 6 x 104, both experiment and 
theory indicate a large a;nount of separation. It seems remarkable that the 
experimentally determined critical Reynolds number agrees so well with that 
predicted by the theory. 

High lift airfoils can also be designed and analyzed with the program. 

One such airfoil designed by Chen (ref. 10) is shown in figure 13.^ Notice 
that the panel method has predicted some oscillations in the velocity distribu- 
tions. (The occurrence of these oscillations is common for the newer airfoil 
designs as opposed to the older NACA airfoils for which the panel method 
predicts smooth velocity distributions.) The objective of this airfoil design 
was to achieve on the upper surface a certain length of constant velocity 
followed by a Stratford pressure recovery (ref. 11). The boundary layer 
development for this airfoil showed early transition due to the oscillations 
in the velocity distributions. These results, of course, were unrealistic and, 
accordingly, a new airfoil, 1220, was designed with the same objective (fig. 14). 
To demonstrate that the oscillations in the velocity distributions for the 
Chen airfoil were not produced by the panel method, the velocity distributions 
from the panel method are included for the new airfoil. 

The boundary layer results for this airfoil a»*e quite interesting. If 
the transition point is just ahead of the pressure recovery, the predicted 
boundary layer remains attached until the closure contribution is reached. 

This occurs at a Reynolds number of 6 x 106 (fig. 15). This demonstrates that 
the method predicts the boundary layer development for an extreme pressure 
recovery quite well. (This had already been tested by M. Schulz.) At a 
Reynolds number of 3 x 106, an intense laminar separation bubble was predicted 
at the beginning of thp nressure recovery. At a Reynolds number of 9 x 10®, 
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transition is predicted further ahead of the pressure recovery and a thicker 
(turbulent) boundary layer arrives at the beginning of the pressure recovery 
than if the flow had remained laminar up to that point. And so, again the 
precise initial conditions for the Stratford pressure recovery are not satis- 
fied and early turbulent boundary layer separation is again the result. This 
must be true for every airfoil having a Stratford pressure recovery derived 
from one initial condition. 


Accordingly, a new airfoil was designed which would not exhibit the 
undesirable characteristics of the previous airfoil. Because maximum lift 
normally occurs, in flight, at lower Reynolds numbers, the objective of the 
new design was to develop a high lift coefficient at a Reynolds number of 
1 X 10® while still maintaining a soft stall for practical reasons. The air- 
foil which resulted is shown in figure 16. The a' distribution for the 
forward portion of the upper surface was chosen such that no sudden movement 
of laminar separation or transition is possible. This feature is demonstrated 
in figure 17 in which all the curves show decreasing with increasing R^ 
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The pressure recovery is concave but not nearly as extreme as the Stratford 
distribution. It should be mentioned that the total amount of pressure 
recovery for this distribution is only slightly less than that of the Stratford 
distribution. Thus, the moderate pressure recovery results in more lift. 


Another feature of this design is that the upper and lower surface 
velocities aheau of the closure contribution were not required to be equal as 
in the case of the Chen airfoil. It has already been demonstrated by the NACA 
6-series airfoils that this condition is not necessary. 


The theoretical results for this new airfoil =»re shown in figure 18. The 
maximum lift coefficients to be achieved by such airfoils are surely above 2. 


CONCLUSIONS 

The present program system of combined potential flow and boundary layer 
theories has been discussed. Applications and comparisons with experiments 
over a very wide range of Reynolds numbers have been shown. The results are 
most satisfactory and open the door to the tailoring of airfoils for specific 
objectives. 
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COMPONENT EQU..LS ZERO AT P AND P 

n n+1 


Figure' 2.- Panel method. 
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WAKE 1: -Pj (KUTTA-CONDITION) 

TRIANGULAR WAKE WITH CONSTANT OPPOSITE VORTICITY 
NORMAL FLOW CONDITION AT Pj AND P 
NO FLOW CONDITION AT THE WAKE ^ 



WAKE 2; ^ -p^ (KUTTA) 

LINEAR VORTICITY AND SOURCE DISTRIBUTION AT THE BASE SUCH THAT NO FLOW 
SINGULARITY EXISTS AT P^ AND P^. FLOW CONDITIONS: INNER TANGENTIAL 
AND INNER. NORMAL VELOCITY EQUALS ZERO AT THE MIDDLE OF THE BASE 


Figure 3.- Blunt trailing edge. 
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Figure 12.- Comparison of experiment and theory for EPPLER 387 airfoil. 






















